PLOS Genetics

Coupling histone gene expression to S phase of the cell cycle is essential for genome duplication and stability. Activation of Cyclin E/Cdk2 at the G1-S transition stimulates high-level expression of histone genes during S phase, but how histone genes are turned off at the end of S phase is not understood. Here we demonstrate that the essential Drosophila gene mute functions to repress inappropriate histone mRNA accumulation outside of S phase by counteracting Cyclin E/Cdk2-dependent phosphorylation of Mxc, which activates histone gene expression. Additionally, Mute plays contrasting roles in histone gene expression during S phase by promoting high levels of H1, H2a and H2b expression but not H3 and H4. Although Mute is present only at replication-dependent histone genes, its loss leads to 801 differentially regulated genes, primarily those involved in muscle related processes in late-stage embryos. Thus, disruptions of histone gene expression control alters the transcriptome resulting in developmental defects.

The continuous renewal of healthy epidermis depends on the finely regulated proliferation of basal keratinocytes and subsequent differentiation as the newly formed cells move upwards through the different layers of the epidermis. Perturbations in keratinocyte differentiation may lead to cornification disorders. We investigated seven dogs of different breeds belonging to four independent families that showed striking multifocal tree bark-like skin lesions. Histopathologically, lesional skin was characterized by pronounced epidermal and infundibular hyperkeratosis with epidermal and sebaceous gland hyperplasia. We therefore tentatively termed the phenotype phloiokeratosis, derived from the Greek word phloios for tree bark and keratosis indicating abnormal keratinization. Whole genome sequencing of DNA from affected dogs revealed four independent variants in the SUV39H1 gene encoding the SUV39H1 histone lysine methyltransferase, an H3K9 methyltransferase, which is involved in epigenetic silencing of chromatin. Phloiokeratosis is inherited as an X-chromosomal semi-dominant trait. Four of the affected dogs in our study were heterozygous females and had lesion patterns reminiscent of Blaschko lines. In two of them, trio analyses experimentally confirmed de novo mutation events in the SUV39H1 gene. Previously, Suv39h1-/- knockout mice had been reported to have normal skin. So far, no human patients with SUV39H1 loss-of-function variants have been reported. The findings in SUV39H1 mutant dogs with phloiokeratosis for the first time link SUV39H1 deficiency to a heritable skin phenotype. Our study highlights the essential role of SUV39H1-mediated epigenetic silencing during normal keratinocyte differentiation and provides a unique model for further investigations. Author SummaryThe integrity of the skin depends on a balanced equilibrium of keratinocyte proliferation, differentiation, and sloughing of terminally differentiated cells into the environment requiring finely regulated changes in the global transcriptome of differentiating keratinocytes. We investigated seven dogs belonging to four different families with a new disorder of cornification characterized by tree bark-like outgrowths of the epidermis. Histopathological examinations confirmed that the outermost layer of the epidermis was thickened in affected dogs. The genetic analysis yielded four different SUV39H1 loss-of-function variants in the affected dogs from the four families. The SUV39H1 gene encodes an enzyme that is involved in the epigenetic silencing of chromatin. The newly characterized inherited skin disease in dogs is the first clinical phenotype that has been linked to SUV39H1 deficiency. Most likely, SUV39H1 deficiency leads to delayed epigenetic silencing and consequently delayed differentiation of keratinocytes. Dogs with this rare skin disease provide an improved understanding of the essential role of SUV39H1 in the epigenetic control of gene expression in skin.

A mild impairment of mitochondrial function activates the hypoxia inducible factor (HIF-1)-mediated hypoxia stress response pathway leading to a HIF-1-dependent increase in lifespan. Lifespan extension resulting from HIF-1 stabilization is dependent on activation of flavin-containing monooxygenase-2 (FMO-2). In this work, we explored the role of fmo-2 in the long lifespan of genetic mitochondrial mutants in C. elegans. We found that fmo-2, but not other fmo genes, are specifically upregulated in the long-lived mitochondrial mutants clk-1, isp-1 and nuo-6. Disruption of fmo-2 through RNA interference or genetic mutation shortens the lifespan of these mitochondrial mutants indicating that fmo-2 is required for lifespan extension in these worms. Moreover, signaling molecules that have been shown to be involved in upregulation of fmo-2 are also required for the long life of clk-1, isp-1 and nuo-6 mutants including HLH-30, NHR-49 and MDT-15. Finally, we examined the effect of multiple lifespan-promoting pathways in clk-1 mutants on the expression of fmo-2. We found that in all cases, genes required for clk-1 longevity are also required for the upregulation of fmo-2 in clk-1 worms. These genes included DAF-16, PMK-1, SKN-1, CEH-23, AAK-2, HIF-1 and ELT-2. Combined, this work advances our understanding of the molecular mechanisms contributing to longevity in the long-lived mitochondrial mutants and identifies FMO-2 as a common downstream effector of multiple pathways that modulate longevity.

Meiotic recombination is a key driver of evolution in sexually reproducing organisms, reshaping genetic diversity by generating novel allelic combinations. The rate of recombination varies substantially across living organisms depending on cis- or trans-acting genetic elements, as seen in many species, including the yeast Saccharomyces cerevisiae. Here, we report on an experimental evolution-based study to better understand the factors shaping this natural variation. Starting with a genetically diverse population of S. cerevisiae, we have carried out recurrent divergent selection on recombination rate using a fluorescence-based sorting approach in four independent lineages. After ten generations, we observed an average response of recombination rate of +28% after positive selection and -24% after negative selection, within the interval used for selection. In the adjacent region, however, we observed a weaker response in the opposite direction, and no response in four other unlinked genomic regions. Whole-genome sequencing of individuals selected for high recombination revealed mixed outcomes in the four independently evolved lineages. All four lineages showed selection for high recombination locally, with particular haplotypes heavily favored and sequence- or structural variation-based heterozygosity selected against within the selection interval. However, only two of the four lineages showed increases in genome-wide recombination rate. Overall, this experimental evolution approach provides original and useful insights into the evolvability of the meiotic recombination rate and the associated genetic determinants.

Translation reinitiation (REI) is one of the most important gene-specific regulatory mechanisms by which eukaryotic cells influence expression of main translons, for example during highly conserved integrated stress response (ISR). In S. cerevisiae, expression of the key stress response gene, GCN4, is controlled by an intricate interplay among four short upstream translons (uTranslons, formerly uORFs), resulting in high or low levels of REI at GCN4 depending on the growth conditions. Under nutrient rich conditions, GCN4 expression is repressed, but upon amino acid starvation, it is derepressed, despite of a general translational shut down. Capitalizing on our screening reporter system, we identified three new factors influencing efficiency of REI after translation of GCN4 uTranslons: Rai1p (an RNA quality control and processing factor), and Ssz1p and Zuo1p (members of the Ribosome Associated Complex [RAC]). Importantly, we showed that depletion of these factors deregulated derepression of Gcn4p synthesis under starvation. Furthermore, we found that similar to RAC, Rai1p associates with 40S subunits and actively translating ribosomes. We also explored interactomes of these proteins. Collectively, we present three previously unknown factors that co-regulate stress response to amino acid starvation in the budding yeast by unique mechanisms.

Stress granules (SGs) are essential subcellular assemblies that enable cells to adapt to environmental stress, and their dynamic assembly and disassembly are critical for maintaining cellular homeostasis and reproductive capacity. However, the function and regulation of stress granules in germline cells remain mysterious. In this study, we characterized the function of GTBP-1, the homolog of human G3BP1, in response to heat stress in Caenorhabditis elegans and elucidated its regulation in stress granule dynamics. gtbp-1 mutants exhibit pronounced temperature-dependent biodirectional reproductive characteristics. At lower temperatures, their brood size is higher than that of wild-type animals, whereas at 25 {degrees}C they are completely sterile. GTBP-1 is diffusely distributed at normal culturing temperatures but forms perinuclear stress granules upon heat shock in the germline. While the NTF2 domain is essential for germline stress granule formation, other domains, including the IDR, RRM and RGG, are required for the recovery phase after heat shock. GTBP-1 stress granules colocalize with the P-bodies, but not with other germ granules. The depletion of P-bodies prohibited the perinuclear GTBP-1 stress granule formation. A combination of forward genetic screening together with RNAi-based candidate screening identified the P body components, the RNA helicase LAF-1, the SUMO protein SMO-1, and the mTOR pathway effector RSKS-1(S6K) as key regulators of GTBP-1 germline stress granule formation. Together, this work revealed that germline stress granules may be subjected to multiple layers of regulations and GTBP-1 may safeguards reproductive homeostasis under temperature stress by coordinating germline stress granule condensation. Author SummaryReproduction is particularly vulnerable to environmental stress, yet germ cells must remain functional to ensure fertility. How germ cells protect themselves under stressful conditions is still not well understood. In this study, we used the nematode Caenorhabditis elegans to examine how germ cells respond to heat stress, focusing on GTBP-1, a conserved protein involved in stress granule formation. We found that GTBP-1 plays contrasting roles depending on temperature: under normal conditions it restrains reproduction, whereas under heat stress it becomes essential for maintaining fertility. When animals are exposed to elevated temperature, GTBP-1 rapidly forms granules around the nuclei of germ cells, these granules assemble at sites occupied by processing bodies, another type of RNA-containing structure, revealing a close spatial relationship between two stress-responsive compartments. We also found that this process is regulated by conserved factors involved in RNA regulation and nutrient-responsive signaling. Together, our findings show how germ cells reorganize RNA-protein assemblies to preserve reproductive capacity under stressful conditions.

The CenIR regulatory system of Clostridioides difficile comprises a predicted transcriptional repressor, CenI, and a predicted membrane metalloprotease, CenR. The physiological role of CenIR and activating signal(s) are not known. CenIR belongs to the BlaIR family of regulators that mediate resistance to {beta}-lactam antibiotics. In canonical BlaIR systems, binding of a {beta}-lactam to the extracellular transpeptidase domain of BlaR triggers proteolysis of BlaI and thus induction of a closely linked {beta}-lactamase gene. However, CenR lacks a {beta}-lactam-binding domain and transposon mutagenesis indicated CenI is essential for viability even when {beta}-lactams are not present. Here we confirmed essentiality of CenIR and determined its regulon contains [~]12 genes, including an exported protein of unknown function (CDR_0474) that is induced about 500-fold and a peptidoglycan hydrolase (Cwp6) that is induced about 7-fold when cells are depleted of CenIR. There are no essential genes or {beta}-lactamases in the regulon. Phenotypic characterization of CenIR-depletion strains revealed slower growth, mild elongation and cell lysis. Deletion of cdr_0474 corrected all three defects, while deletion of cwp6 only rescued the lysis phenotype. It was possible to delete cenIR in either a {Delta}cdr_0474 or {Delta}cwp6 background. We propose that CenIR is essential because its absence leads to lysis due to Cwp6 overproduction. Bioinformatic analyses revealed the predicted extracellular sensing domains in annotated "BlaR" proteins are diverse. Thus, BlaIR systems are not dedicated to defense against {beta}-lactams but probably enable bacteria to adapt to a variety of environmental stimuli. ImportanceMany of the regulatory systems for controlling cell envelope biogenesis and stress responses have yet to be studied. Here we characterize a Clostridioides difficile BlaIR-like regulatory system that we have named CenIR for cell envelope. Unlike canonical BlaIR systems, which bind {beta}-lactams and induce a {beta}-lactamase, CenIR lacks a {beta}-lactam binding domain and is essential for viability even in the absence of antibiotics. We identified the genes in the regulon and found that CenIR is essential because its absence leads to overproduction of the Cwp6 peptidoglycan hydrolase. We also show that most annotated BlaIR-like systems lack a {beta}-lactam-binding domain, from which we infer that these systems have much broader physiological roles than generally appreciated.

LIN-12/Notch signaling regulates C. elegans vulval development via cell fate specifications in the gonad and epidermis. In the somatic gonad LIN-12/Notch activity specifies the anchor cell (AC) versus ventral uterine cell (VU) fates, with VU receiving more signal. The AC secretes epidermal growth factor (EGF) which induces the underlying vulval precursor cells (VPCs) to adopt vulval fates. In the VPCs the secondary vulval fates are specified by LIN-12/Notch activity. We previously reported that the AGEF-1, an Arf GEF homologous to ArfGEF1 and ArfGEF2, the ARF-1 GTPase, and the adaptor protein complex 1 (AP-1) inhibit LET-23/EGF receptor (EGFR) signaling in the VPCs by antagonizing LET-23/EGFR basolateral localization. Here we report that AGEF-1, ARF-1 and AP-1 regulate LIN-12/Notch signaling during somatic gonad and vulval development. The lin-12(n302) partial gain-of-function causes a potent Vulvaless phenotype due to a lack of AC specification. We demonstrate that loss of AGEF-1, ARF-1 or AP-1 restored the AC fate in lin-12(n302) animals, indicating that AGEF-1/ARF-1/AP-1 promotes LIN-12/Notch signaling in the somatic gonad. Interestingly, loss of AGEF-1, ARF-1 or AP-1 also induced ectopic vulval secondary fates in lin-12(n302) animals, indicating that AGEF-1/ARF-1/AP-1 inhibits LIN-12/Notch in the VPCs. Using a LIN-12/Notch biosensor we demonstrate that loss of UNC-101/AP-1 results in decreased signaling in the VU cell and increased signaling in the VPCs that correspond with decreased expression levels of LIN-12/Notch and LAG-1/DSL ligand in the presumptive AC and VU while also causing increased apical localization of LIN-12/Notch in the VPCs. We hypothesize that the differential regulation of LIN-12/Notch signaling could reflect different trafficking pathways in epithelial cells (VPCs) versus non-epithelial cells (AC and VU). Our results indicate that the AGEF-1/ARF-1/AP-1 trafficking pathway maintains the VPC cell fate patterning by limiting both LET-23/EGFR and LIN-12/Notch signaling. Author summaryCell signaling and membrane trafficking are highly interconnected processes whereby membrane trafficking can regulate signal transduction pathways and vice versa. We previously demonstrated that the ARF-1 GTPase, the downstream AP-1 clathrin adaptor and upstream activator AGEF-1 antagonize the membrane trafficking of the Epidermal Growth Factor Receptor (EGFR) and hence signaling during C. elegans vulva induction. Strong loss of the ARF-1 GTPase pathway resulted in ectopic vulval induction. Here we demonstrate that the ARF-1 GTPase pathway differentially regulates Notch signaling to regulate vulva induction. In the somatic gonad it promotes Notch signaling to regulate the specification of the anchor cell which secretes the inductive signal. In the vulva precursor cells, the ARF-1 GTPase pathway antagonizes Notch signaling which cooperates with EGFR signaling to induce the vulval cell fates. We hypothesize that the differential regulation of Notch signaling by the ARF-1 GTPase pathway could be a result of more complex membrane trafficking pathways in polarized epithelial cells (vulva precursors) versus non-epithelial cells in the developing somatic gonad. Thus, the AGEF-1/ARF-1/AP-1 antagonizes both EGFR and Notch signaling in ensuring that only three of the six vulval precursor cells adopt are induced.

Cellular metabolism and gene transcription are closely linked. The conserved transcriptional regulator SIN3 acts as a scaffold for histone deacetylase (HDAC)-containing complexes and is crucial for development, stress resistance, and overall organismal health. SIN3 regulates metabolic gene expression in Drosophila cultured cells, however, an understanding of the extent of its role in coordinating responses to metabolic stress in whole organisms is incomplete. In this study, we explored how SIN3 controls glycolytic gene expression across developmental stages and under genetic and dietary disruption of glycolysis in Drosophila melanogaster. Focusing on four key glycolytic enzymes: phosphofructokinase (Pfk), enolase (Eno), pyruvate kinase (Pyk), and pyruvate dehydrogenase beta (Pdhb), we found that reducing Sin3A levels increases their expression in both larvae and adults, indicating that SIN3 plays a consistent role in balancing metabolic gene transcription. Genetic interaction experiments indicate that Sin3A interacts with Pyk and Eno, regulating transcription in a gene-specific manner. Disrupting glycolysis via genetic or dietary means alters glycolytic gene expression, and SIN3 modulates this response. These findings indicate that SIN3 functions as a metabolic sensor, regulating transcription in response to cellular metabolic stress. Additionally, we demonstrate that reducing Sin3A levels shortens Drosophila lifespan on both low- and high-sucrose diets, emphasizing the importance of SIN3 in longevity. Overall, these results show that SIN3 is a context-dependent regulator of glycolytic gene expression and lifespan in Drosophila, integrating metabolic signals with chromatin-based transcriptional regulation. SummaryTo survive and thrive, organisms must adapt to distinct metabolic inputs. We investigated the response of the conserved transcriptional regulator SIN3 to metabolic stress and its control of glycolytic gene expression in Drosophila melanogaster. By measuring glycolytic gene expression, testing genetic interactions, and assessing lifespan under genetic and dietary perturbations, we found that Sin3A knockdown elevates glycolytic gene expression in a gene-specific manner and decreases longevity. SIN3 also modulates transcriptional responses to disrupted glycolysis and influences lifespan under sucrose stress. These findings identify SIN3 as a context-dependent transcription regulator that links gene expression with organismal metabolic adaptation.

Escherichia coli is highly sensitive to acid and osmotic stress but adapts by modulating the expression of stress responsive genes. Nucleoid-associated proteins (NAPs) play key roles in DNA organization and sensing environmental changes. The histone-like nucleoid structuring protein H-NS is an NAP acting as a global regulator of stress genes. H-NS may alter local chromatin structure to modulate the expression of such genes in response to environmental stress. The H-NS homolog StpA co-regulates several target genes, but its precise role is poorly defined. To investigate the regulatory interplay between these two proteins, we examined transcription, DNA binding and chromatin structure at two regulated operons, hdeAB and proVWX, in E. coli following exposure to acid and salt shock. Our results show that H-NS senses pH and osmotic cues to remodel chromatin and relieve repression, while StpA compensates for H-NS loss, particularly at proVWX, highlighting a coordinated regulatory network.

Changes in regulatory sequences controlling the timing and activity of gene products underlie much of natural phenotypic variation. Yet, what these changes are and how they impact gene expression remain largely unknown. To address this question, we investigated how transcriptional activity and homeostatic responsiveness of orthologous promoters of the metabolic gene TDH3 evolved among Saccharomyces yeast. We found that promoter expression level increased specifically in the S. cerevisiae lineage and that a substantial part of this increase was caused by genetic variants located between the well-characterized, conserved binding sites for two direct transcriptional regulators. These nucleotide changes altered the promoters expression levels while leaving the expression dynamics conserved. Further, the effects of these nucleotide changes were only seen in the presence of a third transcription factor, TYE7p, which is recruited by the other transcription factors through protein-protein interactions. These results suggest that the cis-regulatory changes act through their influence on the collective assembly/activation of the transcription factors, and that changes acting through such a mechanism can allow distinct parts of gene expression, such as expression level and dynamics, to be tuned separately.

The transporting epithelial tissues comprising the Drosophila melanogaster alimentary and renal systems are known to share a common set of enriched genes sometimes referred to as the "epitheliome", reflecting their shared transport functions. Core amongst these genes are the vha genes, which encode subunits of the large Vacuolar-type ATPase (V-ATPase) proton pump complex. However, many of the non-vha components of the epitheliome remain broadly uncharacterised. Here, we explore the role of RNAseK, a gene identified during unbiased epithelial screens in Drosophila whose function within insects is not yet known, though evidence from mammalian systems suggest a role in supporting proton pump activity. We demonstrate computationally that RNAseK is strongly conserved across evolutionary history, and that expression is regulated by the highly epithelial-specific dCLEAR motif. Seeking to understand why epithelial expression is so emphasised, we have assayed the effects of RNAseK knockdown in different epithelia throughout the fly. Across hindgut, midgut, and Malpighian tubules, we note profound defects in gross tissue morphology, transport activity, and fly survival. Mechanistically, our findings that RNAseK may co-localise with the V-ATPase complex, and that V-ATPase inhibition phenocopies RNAseK knockdown, suggest that RNAseK is a critical component of the proton transport axis across Drosophila tissues. Summary statementRNAseK is enriched throughout, and required in, Drosophila melanogaster epithelial tissues. Molecular evidence and evolutionary inferences suggest this is due to a role in the proton transport axis.

Deciphering the genetic architecture of complex quantitative phenotypes remains challenging in quantitative genetics. These traits not only depend of multiple genetic factors but are also established over time and environments. Although quantitative genetics has investigated the genetic determinism of phenotypic plasticity in contrasted environmental conditions, the time related phenotypic plasticity has received less attention. Here we proposed a multivariate Bayesian framework, the Bayesian Varying Coefficient Model, designed for analysing the genetic architecture of the time related phenotypic plasticity by a multilocus approach. We applied the BVCM to time series phenotypes measured at various time scales (daily, monthly, yearly) across a diverse set of biological species. We included in this study: yeast (Saccharomyces cerevisiae), fungi (Fusarium graminearum), eucalyptus (Eucalyptus urophylla x E. grandis), and sweet cherry tree (Prunus avium). The BVCM results were compared with those obtained with a known genome-wide association method carried out time by time. For all species and traits, the BVCM was able to detect the major QTL identified by marker-trait association methods and revealed additional genetic regions of weak effect. It also increased the phenotypic variance explained for most of the phenotypes considered. It revealed dynamic QTLs with transitory, increasing or decreasing effects over time. By considering both the temporal and genetic multivariate structures in a single statistical model, we increased our understanding of the genetic architecture of complex traits notably by reducing the issue of missing heritability. More broadly, this work raises the foundation for extended applications in functional genomics, evolutionary ecology, and crop breeding programs, in which time-related phenotypic plasticity remains crucial for predicting and selecting key quantitative complex traits. Key messageBy capturing the genetic factors influencing the time related phenotypic plasticity, our approach contributes to a deeper understanding of the dynamic nature of genotype-phenotype relationships.

Translation elongation and efficiency are modulated by the genetic code. In the yeast Saccharomyces cerevisiae, 17 inhibitory codon pairs, distinguished by requirements for wobble decoding and distinct codon order, result in reduced translation efficiency and slow translation. Nine of these inhibitory pairs are functionally important as they are disproportionately strongly conserved within the orthologous genes in Saccharomyces sensu stricto. For three pairs, including CGA-CGA, inhibition is triggered by ribosome collisions and known quality control responses, but the mechanisms by which other pairs cause inhibition is unknown. Here, we examined of the molecular basis of inhibition by the slowly translated, highly conserved Leu Pro CUC-CCG codon pair yielding four findings. First, inhibition is mediated by tRNALeu(UAG), which decodes CUC by a U[bullet]C wobble interaction and effectively competes with the nonessential W[bullet]C base pairing tRNALeu(GAG). Second, despite nearly universal conservation of U33 in tRNAs, the C33 alteration in tRNALeu(GAG) does not significantly impair its function. Third, inhibition is likely due to ribosome collisions as many suppressors have mutations predicted to reduce ribosome concentration, including mutations in large ribosomal subunit proteins, RNA polymerase I, and ribosome assembly factors. Furthermore, local reduction in ribosome concentration suppresses inhibition. Fourth, we find a link between the metabolic state and CUC-CCG inhibition, as we find six suppressor mutations in SCH9, a downstream effector of TORC1 that mediates ribosome production. As Sch9 is inactive during starvation, causing reduced ribosome concentration, one biological function of inhibitory pairs may be to mediate a change in relative expression during starvation conditions.

The actin polymerization machinery, comprising the ARP2/3 complex and its activators, the WASP family proteins, has been implicated in regulating a broad spectrum of nuclear processes, such as transcriptional regulation and nuclear organization. Here, using clustered protocadherin (cPcdh) and {beta}-globin genes as model systems, we showed that WAVE2, a member of the WASP family, regulates chromatin organization by maintaining heterochromatin dynamics. Specifically, by CRISPR DNA-fragment editing, in conjunction with integrated analyses of ChIP-seq, MeDIP-seq, ATAC-seq, 4C-seq, and RNA-seq, we showed that deposition of H3K9me3, a key heterochromatin mark, is significantly decreased at the cPcdh locus upon WAVE2 deletion, concurrent with aberrant accumulation of CTCF/cohesin complex at promoter regions and spatial reorganization of chromatin architecture around nucleolus. In addition, REST/NRSF exerts a similar heterochromatindependent effect on the cPcdh locus. Finally, genetic and genomic data showed that WAVE2 regulates {beta}-globin gene expression by maintaining heterochromatin status. Together our data suggested that WAVE2 and REST/NRSF regulate clustered gene expression in a heterochromatin-dependent manner.

Text abstractsLipid homeostasis is essential for organismal physiology, and its disruption contributes to metabolic disorders. Using an unbiased genetic modifier screen in Drosophila, we identified GAR1, a core component of the box H/ACA small nucleolar ribonucleoprotein complex, as a pivotal regulator of systemic lipid storage. We show that the H/ACA snoRNP complex is essential for maintaining lipid droplet morphology in adipose tissue and preventing ectopic fat accumulation. Moreover, null mutants of Gar1 or Dkc1 exhibit severe developmental defects, including reduced body size and larval lethality. RNA-seq analysis revealed that Gar1 dysfunction triggered widespread alternative splicing defects, specifically targeting key transcripts within the insulin signaling cascade, including chico, Pi3K92E, sgg, and Lip4. Furthermore, knockdown of Gar1 impaired insulin signaling, as evidenced by the reduced membrane localization of the tGPH fluorescence. Genetic epistasis further positions GAR1 upstream of the lin-28/foxo axis, as knocking down lin-28 or foxo fully rescues the lipometabolic defects in GAR1-deficient animals. These findings reveal a previously unrecognized link between the snoRNP machinery and metabolic process, establishing the box H/ACA complex as an important coordinator that integrates RNA processing with insulin-mediated nutrient sensing to ensure developmental and lipid homeostasis. Article summaryLipid metabolism is tightly controlled by multiple factors. To find new regulators, the authors performed a genetic screen and identified a small nucleolar protein GAR1 participate in fat storage and larval development. They demonstrated a critical role of box H/ACA snoRNP complex in modulating alternative splicing and balancing insulin cascade. Blocking two insulin-related genes reversed the lipid defects caused by Gar1 loss. These findings revealed the box H/ACA complex integrates RNA processing with insulin-mediated nutrient sensing to ensure developmental and lipid homeostasis, offering a perspective for understanding the metabolic regulation network.

Detecting signatures of positive selection in genomes is a common application of population genetics and one of the most influential models for this task is the hard selective sweep where a de novo mutation rapidly fixes. Many statistics have been developed to detect hard sweeps, often attempting to summarize signatures left behind in the site frequency, spectrum, linkage disequilibrium, and haplotype frequency. However, potentially undiscovered signals could still exist. We attempted to test whether any undiscovered signatures of hard sweeps exist by comparing machine learning models, which can learn signatures from raw data without any prior knowledge, to known summary statistics for inferring the time to fixation (tf) of a hard sweep in a background of variable sweep ages (ta). Across approximately 200,000 simulations encompassing 5 different demographic scenarios of single panmictic populations, machine learning models trained directly on raw genotype data failed to better predict tf than methods based purely on common summary statistics. This suggests few undiscovered signals remain in single timepoint, single population genotype data that can better disentangle tf and ta of hard sweeps.

Gene amplification is a potent driver of evolution and is thought to contribute to genetic diseases, including cancer. The yeast Saccharomyces cerevisiae is a powerful organism for understanding amplification mechanisms. When yeast is grown long term in sulfate-limiting chemostats, amplification of the gene that encodes the primary sulfate transporter, SUL1, is a common outcome. Here we describe a form of SUL1 amplification in which multiple copies of the right terminal region of chromosome II are appended in tandem to a native telomere. We find this form of amplicon when we delete the origin of replication next to SUL1 or delete a variety of genes involved in DNA metabolism. It is the only form of amplification found in a yku70{Delta} mutant suggesting that unprotected telomeres are involved. We propose that these terminal addition events occur when the unprotected 3 G1-3T telomeric sequence invades a short ([~]7 bp) internal telomere sequence (ITS) to begin a form of microhomology-mediated break-induced replication (mmBIR) that has been documented in type-I survivors of telomerase mutants. In addition to amplification of the right end of chromosome II we also find that telomeres containing the sub-telomeric repeat Y experience similar tandem amplification events and show that their formation is reduced in a pol32{Delta} mutant, a gene required for mmBIR. Within individual amplicons the ITSs and Ys are nearly identical, suggesting that the multiple copies of the amplified region are generated in a single mmBIR event that we describe as pseudo-rolling circle mmBIR. A similar amplification event at the P-telomere of human chromosome 18 has four copies of a [~]54 kb region separated by ITSs of nearly identical size. This finding suggests that these additional copies of the terminal fragment of human chromosome 18 arose by the same pseudo-rolling circle mechanism, perhaps during a period of telomeric stress. AUTHOR SUMMARYThe human genome is peppered with duplicates (or higher numbers) of segments that are located at sites both nearby and distant from the original, ancestral segments. These Copy Number Variants, or CNVs, appear to be highly variable among different individuals and are being examined with great interest as potential loci associated with genetic disease. Experimentally determining how these CNVs arise and become distributed across the genome is nearly impossible using humans. We are using budding yeast as the model organism to explore mechanisms of gene amplification. In this work we show that by destabilizing the ends of yeast chromosomes (telomeres) or by interfering with genes involved in the replication, repair, or recombination of DNA results in a specific form of segmental copy number increase that is initiated at telomeres. We propose that a telomere invades an internal chromosome site and sets up a pseudo-circular template for conservative DNA replication. The outcome is a chromosome with multiple, identical copies of a chromosome end arranged in tandem. We believe that it is also a major mechanism used by cells to repair telomeres that have become eroded during aging.

Eukaryotes typically maintain telomere length within a defined range. While short telomeres are known to activate DNA damage responses and limit cell proliferation, long telomeres are associated with extended proliferative capacity. The broader cellular consequences of long telomeres are comparatively less well understood. In budding yeast Saccharomyces cerevisiae, long telomeres have been shown to influence gene expression at specific loci, but whether long telomeres affect transcription genome-wide has not been reported. Here, we analysed transcriptomes in a lineage that inherited long telomeres (originally due to a rif2{Delta} mutation). Transcriptomes were assessed over two rounds of mitosis and meiosis in the absence of the rif2{Delta} mutation. We show that strains with long telomeres exhibit a distinct gene expression profile, including upregulation of membrane transporters and downregulation of a smaller subset of genes. Both up- and down-regulated genes were distributed across the genome, arguing against a purely telomere-proximal effect on gene expression. Affected genes were enriched for Rap1 binding sites, consistent with a model in which long telomeres sequester telomere-associated transcriptional regulators, such as Rap1, and thereby affect gene expression at non-telomeric binding sites for these regulators. Accordingly, the magnitude of transcriptional changes was greatest in strains with the longest telomeres. Together, our findings demonstrate that long telomeres induce a genome-wide transcriptional response that can accompany inherited long telomeres across generations. Similar effects of long telomeres are likely to occur in other eukaryotes, including humans, where long telomeres are associated with disease. Article summaryTelomeres protect chromosome ends, and their length is tightly regulated. While short telomeres are known to be harmful, the effects of long telomeres are less well understood. Using budding yeast, we show that inherited long telomeres alter the expression of dozens of genes across the genome, particularly membrane transporters. These changes are consistent with a model in which long telomeres sequester regulatory proteins away from other loci. Our findings may have broader implications in more complex organisms, including humans.

Genetic regulation of immune cell composition plays a crucial role in the etiology of complex diseases, yet remains poorly understood. We propose a unified analytical framework that integrates genome-wide association studies (GWAS) of cell type proportions with cell-type-wide association studies (cWAS) to systematically characterize both the genetic regulation of immune cell composition and its downstream effects on disease risk. Using single-cell RNA sequencing data from the OneK1K cohort, we conducted a GWAS of immune cell-type proportions with a depth-weighted quasi-binomial model designed for bounded, overdispersed traits. We identified 47 genome-wide significant loci influencing eight fine-labeled immune cell subtypes. Leveraging these identified genetic effects, we further imputed genetically regulated proportions (GRPs) using polygenic risk score (PRS)-based imputation and assessed their associations with complex diseases through cWAS. We identified five significant cell type-disease associations, including two with type 1 diabetes, two with Crohns disease, and one with ulcerative colitis. Together, our results demonstrate that cell type proportions observed in scRNA-seq can reveal regulatory loci and offer insights into how genetic variations regulate immune cell type proportions to affect disease risk. Although we focused on immune single-cell data, our framework is applicable to other tissues or cellular compositions as scRNA-seq datasets expand. Author SummaryGenome-wide association studies (GWASs) have uncovered many disease-associated signals, yet most lie in noncoding regions and are difficult to interpret. Mapping GWAS signals to the relevant cell types is therefore important for better understanding the biological mechanisms that drive disease. A major challenge is that observed gene expression and measured cell-type proportions can be influenced by environmental factors and disease status. In contrast, genotypes are less affected by these factors, making them more reliable for interpreting factors of diseases. Moreover, the cell-type proportions are bounded and often skewed, so standard GWAS models that rely on Gaussian assumptions may lose power. To address this, we developed a quasi-binomial approach that better matches the data and improves discovery while controlling false positives. In real data, our method identified more genetic loci associated with cell-type proportions than a traditional linear model. To further investigate how genetic variation regulates immune cell composition to influence disease risk, we integrated our results with disease GWAS summary statistics to identify immune cell types that may contribute to disease susceptibility. Together, our results link disease-associated GWAS signals to specific immune cell types and provide insights into the cellular mechanisms that may underlie these diseases.